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Mean square matrix elements of the time reversal invariance violating (TRIV) 
interaction between the compound nuclear states are calculated within the 
statistical model, using the explicit form of the TRIV interaction via the p- 



meson exchange. From the comparison of the calculated values with the data 
known for p + 27 Al <^ 24 Mg + a reaction, and for 7-correlation measurements 
I in 113 Cd(n,^/) llA Cd process, the bounds on the TRIV constant are obtained 

g p < 1.8 x 10~ 2 and g p < 1.1 x 10~ 2 . The sensitivity of the recently proposed 
detailed balance test experiments on isolated resonances in 32 S to the value 
of g p is shown to be as high as to reach values g p ~ 10~ 4 . 
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The search for the effects of direct violation of the Time Reversal Invariance (TRI) 
and the related theoretical studies still remain to be a very actual direction of nuclear 
physics lfrj-1^1, fT3^p3|. Different kinds of experiments (nuclear, atomic and molecular) are 



competing in this field to obtain the tighter bounds on the TRI violating interaction. 

It has been realized 0-0] that the TRIV (P even) nucleon-nucleon interaction, Vj , can 
be well approximated by the p-meson exchange, at least at the tree level. (The contributions 
of the lighter mesons are forbidden by the Simonius theorem f5|, while the effects due 
to heavier mesons should be suppressed by the core repulsion.) This interaction can be 
constructed uniquely [§], |J, @ and it is characterized by the only dimensionless constant 
g p that measures the ratio of the TRI violating pNN amplitude to the normal TRI preserving 
one, g p §. Recently, the parameter g p has been adopted as the natural universal measure 
of the TRI violation in the nuclear processes, while the p-meson-mediated TRIV interaction 
is used in the microscopical calculations of the TRIV observables. 

Haxton, Horing and Musolf |J analyzed the limits on the value of g p resulting from 
the existing experimental information obtained in atomic and nuclear experiments. Engel, 
Gould and Hnizdo ]7[] used the interaction V? to obtain the effective TRI violating 
nuclear optical potential. That enabled them to calculate microscopically the "fivefold" 
TRIV correlation (s- 1 x p)(I ■ p) and, by comparison of the result to the experimental value 
for the neutron scattering on the aligned 165 Ho target f21 ], to extract the limits on the value 
of the constant g p . 

However, the experimental limits on the TRIV interaction scale obtained so far, are 
still several order of magnitude greater than the values one can expect from the theoretical 



considerations |10| , 0. Recently, the work |TT| reported on the results of the TRIV scale 
obtained in from consideration of the effective Lagrangians. According to the results of 
Ref. ||IT[, no TRIV effects could be expected above the level g p ~ 1CT 6 (even the stringent 



estimates were obtained in ]T0[), and this makes a serious challenge to the contemporary 
experimental techniques. 

The aim of this work is to show that the limits on the constant g p that can be extracted 



from the statistical nuclear reaction of Ref. pO] and the 7-ray correlation experiments p3 



arc 



close, and they can compete with the constraints obtained from the atomic dipole moment 
data (g p ^$ 9.3 ■ 10~ 3 ) ||, and to show that in principle much stricter bounds on the value 
of constant g p ^ 10 -4 can be reached in the experiments with compound nucleus 32 S", which 
were proposed by Drake, Mitchell et al. in Refs. ||, 0. 

In the compound nuclear reactions, the symmetry breaking effects are usually enhanced 
Hi) [Hi) [0)H1]; e -g-> m the case of the spatial parity violation |24H28f. The TRIV observables 
in compound states are known to be expressed through the mean square matrix element 
of the TRIV interaction P Jl^JT5| -[T7|. Using the explicit form of the TRIV interaction 
Hamiltonian via the p-meson exchange [0], ||, |fLl|| , we make the microscopic calculation of 
the means square matrix element (MSME) of the Vj RIV between the compound states using 
the statistical method developed in p9| , and express MSME through the constant g p . The 
conclusions of the work are drawn from the comparison of the results obtained for particular 
compound nuclei with the data. 

We use the TRIV nucleon interaction via the p-meson exchange (see Refs. [0], ||), that 
reads 



g 2 fji v e -m p \r 1 -r 2 \ 

Vtriv = gp-^—nin x T 2 } z -, ^(1 + m p \r l - r 2 \)(a 1 - a 2 ) x (f x - r 2 ) ■ (p x - p 2 ) (1) 

87cm z ri — r 2 \ 6 



where g p is the dimensionless strength parameter that measures TRIV pNN coupling in 
terms of ordinary TRI conserving coupling constant g p = 2.79, m p and m are the p-meson 
and the nucleon masses respectively, \x = 3.7 is the isovector nucleon magnetic moment in 
units of the nuclear magneton, fj.,2 and p\ y2 stand for the coordinates and momenta of the 
nucleons 1 and 2. 2 and r 12 are the spin and isospin Pauli matrices. The sign x means 
the external vector product. 

Due to the factor {j\ x r 2 } z , only the exchange terms of the interaction (|1|) are operative. 
It is therefore convenient to rewrite the interaction in the form V^ RIV , so that only the direct 
terms should be taken into account 
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v$biv = ^E^cn^, (2) 

ab,cd 

Using the Fiertz transformation for the spin and isospin matrices, 

{ti x f 2 } 3 -> z(fi - f 2 ) 3 , (cti - <r 2 ) 3 -> -^{cti x ct 2 } 3 
the matrix elements Vtriv can be found in the form 

2 , 



fr/ab,cd — 9p^v f , , 



2ia 



1/3? 



m p e 



[crz a 



f-afa c Q d )Mn)Mr2), (3) 



47r|ri — r 2 | 

where Vi = ipi, dl and d2 mean integration and summation over the spatial and the spin 
variables, [ , ] means the commutator and the summation over the double indices is 
assumed. Like the ordinary (symmetry preserving) strong interaction, the proton-neutron 
interaction ([!]), (|3]) has P-even selection rules, with the essential difference that the matrix 
elements of Vtriv are purely imaginary, and therefore, its diagonal matrix elements vanish. 

1. The method of calculation of mean squared weak matrix element is based on the 
equivalence theorem of canonical and microcanonical ensembles for a system with a large 
number of degrees of freedom. It has been developed previously for the microscopic calcu- 
lation of the parity violating matrix elements between the compound states f29fl . (Later, it 
has been applied to the calculation of the P and T noninvariant matrix elements p0[). We 
recall that the wave function of any compound state with angular momentum J and parity 



7r, | J 71 ), can be expressed as the superposition of the simple components | J 71 " > [[31 



\r) =Y J C a \Jl>, \J%>= (a + bc + de + ...)j,\G.s. > (4) 

a 

Here, \J* > are the multi-particle excitations over the shell-model ground state \G.s. >. 
By (...)jt, we mean the coupling of nucleon creators a + and destructors a to total angular 
momentum J and parity n. The energies E a of the "principal" components \J£ > which 
dominates the normalization of (^), must be within the interval [S — -^7p,£ + ^-rp] 



24[ , where 8 is the energy of the compound state and T spr is the spreading width (typically, 



T spr ~ 1..3MeV, Refs. [0). These components (which contain several excited nucleons) can 
be composed by excitations of protons and neutrons only inside incomplete valence shells. 
Mean squared values of the coefficients C(E a ) can be described by the formulae (see e.g., 
Ref. TO 



C*(E a ) = J7 A(T spr ,S-E a ) } A/-=^p, (5) 

^ij^spri £ E aj 



r 2 /4 

sprl 



(£ — E a ) 2 + r 2 pr ,/4' 

where E a is the energy of an arbitrary many-particle configuration), d is the averaged energy 
distance between the resonances, and J\f ^>1 is the number of principal components. The 
Breit-Wigner-type factor A describes quenching of weights for states distant in energy. It is 
a "spread" 5- function normalized as to be of order unity for \£ — E a \ < r spr /2, and obeying 
the limit A(r spr , £ — E a ) — > nr ° pr 5(£ — E a ) for r spr — > 0. Thus the average with weight C 2 
is equivalent to the microcanonical average. 

The TRIV interaction (Eq.fll])) is a two-body operator (|2|), and it has nonzero matrix 
elements between the single-particle states of the valence shells. Its matrix elements between 
the compound states (Q) can be expressed in terms of the single-particle matrix elements 
Vtriv by using of the statistical properties of the components of the compound states l . 



To calculate the mean squared value of the TRIV matrix element |V^ HJl/ | = v\ Rlv one can 
write it in the form 



VTRIV = I ( J* I Vf RIV \J*')(J*' | V£ mv I J") 



1 The case of the parity-odd interaction considered in [29 is a bit more involved. The most 



important contribution to the parity odd MSME is given by the single-particle symmetry breaking 
potential which acts in combination with the symmetry preserving residual strong interaction to 
dominate the MSME. In the present T-odd case, the corresponding single-particle potential that 
can be obtained from Eqs.(||j2|,^), does not contain coherent contribution from the core nucleons 
0. Therefore, the T-odd MSME is dominated by the matrix elements of ([!]), (|3]). 
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= E \C' a C' p {J-\VS mv \Jl >< J$\Vf RIV \J*)\ (6) 

a/3 

using the expansion of the compound state | J ,7T ) in terms of their multiparticle components 
(|j). Here | J /7r ) and \J n ) are the different compound states with the same quantum numbers. 
The coefficients C a are statistically independent [[H]], C a Cp = C^5 a /3 = 5 a pjjA{Y spn £ — E a ), 
[Eq.(|5])]. Therefore, in (||) we obtain 

VTRIV = E JJ'^C^spr, £ — E a ) x 



x\(J*\v£ RIV \JZ >< JZ\V£ RIV \J*)\, (7) 

In the second quantization representation (0), summation over the intermediate states a 
in the last equation is equivalent to summation over different components of the operator 



Vtriv- m the expression for the remaining average, (J n \ Vj, RIV V^ RIV \ J n ), in (0), the weights 



CC will appear in the same way as above. To calculate this average over the resonances 
\J n ), we now use, instead of the microcanonical ensemble @, the equivalent canonical 
(Gibbsean) ensemble. (As we know from the statistical mechanics, the latter can be used 
as an equivalent counterpart of the microcanonical ensemble, provided M 3> 1.) Thereby, 
the average expectation value in ([?D can be evaluated using the standard contractor rules 



for the Fermi operators, given by Wick theorem ( J n \a + b\ J w ) = (a + b)x = 5 a b^ a (T), for v a (T) 
being the finite temperature Fermi occupation probabilities of the single-particle levels, 
f a (T) = {exp[(e a — A)/T] + l} -1 . The parameters of the effective canonical ensemble, T (the 
temperature), and X p , and A„ (proton and neutron chemical potentials) are to be determined 
from the "consistency " equations E = J^a^a^a, Z = J2 p v p , and iV = J^n^n for the total 
energy E, nuclear charge Z, and neutron number iV correspondingly. Making use of the 
Wick theorem reduces (||) to the result (see fl29f): 



vtriv — \J\Vtri\ 



d 



£ u a (T)[l - v b {T))v c {T)[l - u d (T)} | V TRI f cd | 2 A(T spr ,e a -e b + e c -e d 



\ TlT 



s P r \abcd 



Here, V TRI y ,cd is the direct matrix element of the TRIV Hamiltonian in the form (^), defined 
according to Eq.(§). The factor A reflects the required energy balance in the transitions due 



to V^ RIV and ensures that the components of both the states | J 71 ) and | J /7r ) are in the same 
energy domain £ ± Y spr /2. 

The numerical calculations of the vtriv have been performed with the use of single- 
particle basis of states obtained by numerical calculations of the eigenvalue problem for 
the Woods-Saxon potential with spin-orbital interaction in the form U(r) = —Uof(r) + 
Ui s {al){h/(m^c)f\f r + U c with /(r) = (1 + e (r-R)/*y\ w h ere f is the orbital angular 
momentum, U c means Coulomb correction for protons, U c = 3Ze 2 / (2R)(l—r 2 / (3R 2 )),r < R 
and U c = Ze 2 /r,r > R, for R, a, and r being the nuclear radius, diffusity parameter and 
radial variable correspondingly. The parameter values were used in accordance with Bohr- 
Mottelson prescription [[Hp. As the interaction via the exchange by the massive p-meson 
is affected by the repulsion at the core scale distances, the short-range correlations were 
accounted for within the standard Jastrow method (see, e.g., p2| ). 

2. There were two experiments with reaction p + 24 Mg <(=> 27 Al + a to study TRIV via 



detailed balance test [[HJ, PD| . They give upper limit on the value of the parameter a that 
is the ratio of the r.m.s. of the TRIV matrix elements between the compound states to the 
regular TRI-preserving strong interaction Vg: 

vtriv , n s 
« = 7=^1/2- ( 9 ) 

T/2 ' 

The revision of the data of the improved experiment [20| by Harney, Hiipper and Richter 



1 22|| gave the following limit: 

a < 2.6 ■ 10" 4 (80 % conf.lev.). (10) 

To obtain the theoretical value of the quantity a in terms of the universal constant g p , 
we need also the r.m.s. of the residual strong interaction matrix elements between the 
multiparticle components of the compound states. We have employed the most widely used 
Landau-Migdal particle-hole interaction which rises to Landau Fermi liquid theory. For the 
case of a nucleus it was justified in the Theory of Finite Fermi Systems P5|P3| . In explicit 
form, the Landau-Migdal interaction reads 
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V s = C5(n - r 2 )[/ + f(m) + g(w 2 ) + g'itt)^)], (11) 



where C = 300 MeV fm 3 is the universal Migdal constant [f£| , /, /', (7, g' are the density 
dependent strengths / = / in - (/«,,- fm){p{r)- p(0))/ p(0), p{r) = p(0)/(r). The interaction 
(PI) has been successfully used by many authors (see Refs. pi| ) to describe a great amount of 
various properties of nuclei at quantitative level. The values of the parameters [f ex = —1.95, 
f m = -0.075 f ex = 0.05 f' m = 0.675, g in = g ex = 0.575, and g' m = g' ex = 0.725] were 
chosen to fit a variety of experimental data. It is therefore reasonable to believe that the 
parametrization ( |TT|) describes strength and the main properties of the residual interaction 
correctly. 

The value of Vf can be calculated via the same formula as vtriv © with the 
substitution the direct matrix elements of (JTTj) into (H) instead of V^ RI y' cd , (||). In this case, 



we need in fact to calculate the MSME of V between the components of the compound 
states {{J 71 >} (|4]) rather than between the compound states {|</ 7r )}. However, we can 
use the same formula (^|), as the sum of the squared matrix elements taken in the energy 
interval T spr is independent of reference basis, while the same thermo dynamical expression 
for the occupation numbers v can be used [[!!]] in averaging over the unmixed components. 



In contrast to the TRIV interaction, the normal strong interaction that is given by the 
direct terms of (|TTJ) contain the the diagonal matrix elements (J' K \Vs\J' K ) ■ Those have to be 
excluded from the sum (H). 

In the expression for a (§), all the coefficients in front of the square root of the sum in @ 
are canceled. Using the value of temperature T = 2.3MeV to reach the excitation energy 



E exc ~ 16..20MeV of the compound nucleus 28 Si W (T spr ~ lMeV |20|), we have obtained 



the value of a, and, consequently, the upper limit on the TRIV constant g p , as follows 

a = 1.44 -10-%, g p < 1.8 • 10~ 2 . (12) 

This limit is close to the one obtained in Ref. 0] from the atomic dipole moment data 
(g p < 0.93 • 10~ 2 ). It should be borne in mind that twelve years passed after the experiment 



of Ref. 20 . 



3. We consider now the result of the novel experiments |2l| where the forward-backward 



asymmetry in individual 7-ray transitions resulting from unpolarized neutron capture in 
p-resonances by the target ll3 Cd has been measured. (The possibility of measurement of 
T-odd correlations in the experiments of such type was discussed earlier in Ref. |3S[| .) The 
value of TRIV mixing coefficient (3 has been extracted for the compound nucleus U4 Cd 
(neutron p-resonance at energy E p \ = 7.04eV) to be 



P 



VTRIV 



(3 EXP < 0.08 (13) 



Assuming the mixing with the closest p-resonance (£ p i =21.8 eV H37||), we obtain, for the 
value of TRIV matrix element vtriv ~ 1.18eV. The present calculation of the root MSME 
vtriv f° r the low energy neutron resonances in U4 Cd (the temperature T = l.OAMeV, 
d = 27eV, and T spr = 3MeV) gives the following quantity and corresponding bound 

vtriv = 9 P ■ 0.11 • 10 3 eV, g p < 1.1 ■ 10" 2 , (14) 



This constraint is tighter than in the case of 28 Si ([12]). One should note that the bound 
on g p extracted from /3, depends strongly on the density of the p-resonances. E.g., suppose 
the same upper limit on (3 ([TJj) to be obtained on one of the resonances at £ p = 98.5 eV 
and S p = 102.3 eV | 37fl . This would result in the lower bound on the TRIV constant: 



g p < (2 - 3) x 10- 3 . 

To judge on the reliability of the statistical method used here, it is interesting to see the 
results the method gives for the P-odd, TRI preserving weak interaction for the compound 
systems similar to those considered here, where the experimental data are known. (Earlier, 
the method produced the result consistent with experimental data for 233 Th compound 
nucleus P5f| ). For 28 Si and 32 S (discussed below), the "closest case" is the compound nucleus 
36 CZ, where the experimental value of the P violating MSME vpy = 60(±20)mel / 0] was 
obtained in the low energy resonances. The theoretical result obtained by the present method 
(as described in [J29|), is 26meV. For cadmium, there are data in low energies for the same 
isotope, U4 Cd vpy = 1.4(±0.5)meV. The corresponding theoretical value is found to 
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be 1.2meV, for the same values of d, T spr and T as above. Therefore, the method seems to 
give reliable results. 

4. Now, let us estimate the scale of the constant g p which can be available for in- 
vestigation in the experiments with with formation of the compound nucleus 32 S in the 
energy region of separated resonances, proposed in || (detailed balance test in reaction 
p+ 31 P <^ 28 Si+a). In this case, the observable T-odd quantity is the asymmetry parameter 
A that measures the degree of the detailed balance violation [|8|, A = 2 ffi ''°7' T °' p . Here, o v a 
and the properly weighted differential cross sections for direct and inverse processes, 

see ||. As was shown in Refs. ||, 0, the quantity A is proportional to the average value 
of the matrix element of the T-odd interaction vtriv 

A = T ■ vtriv, 

where the factor T depends in a complicated way on the angles and resonance parameters. 
Its absolute value \T\ = I — - — I determines the sensitivity of the experiments to the T-odd 

I I I v TRIV I J r 

interaction. The analyis of Refs. || , || has shown that high sensitivity can be produced by 
an appropriate choice of the conditions of measurement (energies, angles etc). 

Basing on the the present results for the vtriv [Eq.@] and on the values of \F\ obtained 
in Refs. ||, |§, we can make the estimation for the upper bound on the TRIV constant g p 
that could be reached in such experiments. The values of the resonance energy spacing and 
the spreading width were take from Ref. f9|: d — 7.2 ■ 10 4 eV and r spr = lMeV. The 
value of temperature T = 2.1MeV was used in accordance to the consistency condition 
for excitation energy (see above) to reach the excitation energy E exc ~ 10..12MeV of the 
compound nucleus 0. For the 'best sensitivity' in range | V J^ IV I ~ (1-5 — 15) • 10 -3 0, 0, 
the results are summarized as follows: 
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Level of the Present result Bound on g p 

vtriv expected 

measurable effect @ for the v TRIV as resulted 

g p < 1.4 • (1CT 3 - 

A ~ 1(T 2 6.7 - 0.67 eV gA.8 • 10 3 eV 

lcr 4 ) 

Therefore, it can be concluded that the experiments for detail balance tests of TRI with 
32 S'-compound nucleus proposed in Ref. 0] can produce not only the new bound on the 
TRIV constant g p , but lower the upper limit for g p down to 10~ 3 — 10~ 4 . This is already 
close to the most "pessimistic" expectations for the value of constant g p that follow from 



the considerations by Engel and collaborators (11 



In summary, the mean square matrix elements of the TRIV interacting between the 
highly excited nuclear states have been calculated within the statistical approach, starting 
with the microscopic TRIV nucleon interaction via p-meson exchange. The r.m.s. values 
of the TRIV matrix elements are expressed through the p- mediated interaction constant. 
Particular calculations have been made for the 28 Si, 11A Cd and 32 S compound nuclei. The 
results, being compared with the experimental data from the compound nuclear reactions, 
allow one to obtain the limits on the TRVI constant g p < 1.8 x 10~ 2 and g p < 1.1 x 10 -2 , for 
the detailed balance test experiment in the statistical reaction in the region of overlapping 



resonances [£(]] an d in the 7-correlations in low-energy neutron resonances ||23|1 , respectively. 
The sensitivity of the experiments of the detailed balance tests in the p + 31 P -v^ 28 Si + a 
reaction, proposed in ||, to the value of TRIV p-meson constant is found high enough 
to study the new range of values g p ~ 10 -4 . 

The conclusion can be drawn that the nuclear experiments with compound states seem 
rather prospective and are able to provide important information on the TRIV scale. Of 
course, the TRIV matrix element between the compound states is a fluctuating quantity, 
and hence the bounds obtained from vtriv are approximative. A more definite information 
can be obtained from the measurements on several resonances and/or in different compound 
systems. 
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One of the interesting issues related to the TRIV interaction via the p-meson exchange 
(|I|) is that explicit microscopic form of the TRIV Hamiltonian opens the possibility to 
study correlations between TRIV matrix elements and the matrix elements of the normal 
TRI preserving strong interaction. (Usually, the TRI preserving and TRI violating matrix 
elements have been treated as statistically independent random quantities, i.e., the members 
of the Gaussian unitary ensemble.) Apparently, these correlations should be depressed 
for the same reason as the correlation between P-odd and P,T-odd matrix elements |30| . 
However, this interesting question deserves studies. 

The author is grateful to V.V.Flambaum for his interest to the work and comments, 
V.F. Dmitriev and V.B. Telitsin for kindly providing the author with the codes to solve eigen- 
value problem in Saxon- Woods potential, and to V.A.Kuznetsov for reading the manuscript. 
The work has been supported by ARC grant. 
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